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Abstract. In this paper, a family of random Jacobi matrices, 
with off-diagonal terms that exhibit power-law growth, is stud- 
ied. Since the growth of the randomness is slower than that of 
these terms, it is possible to use methods applied in the study of 
Schrodinger operators with random decaying potentials. A partic- 
ular result of the analysis is the existence of operators with arbi- 
trarily fast transport whose spectral measure is zero dimensional. 
The results are applied to the infinite Dumitriu-Edelman model [6] 
and its spectral properties are analyzed. 



1. Introduction 

For a sequence of positive numbers, {a(n)}™ =1 , and a sequence, 
{b(n)y^ =1 , of real numbers, let J ({a(n)}^ =1 , {b(n)}™ =l ) denote the Ja- 
cobi matrix with off-diagonal elements given by {a,{ri)}™ =1 and diagonal 
elements given by {b(n)}'^ =1 on the diagonal. That is: 



J(Un)}~ =1 ,{b(n)Y f 



( 6(1) a(l) 

a(l) 6(2) a(2) 

a(2) 6(3) a(3) 

V i •• •• ■• 



/ 



(1.1) 



For rji e (0, 1) and A x > 0, let J\ uril be the Jacobi matrix whose 
parameters are a\ 1M {n) = \in Vl and b\ ltVl (n) = 0. Since rji < 1, J m 
is a self-adjoint operator on £ 2 (N) [1]. For any such operator we may 
define the spectral measure, /i, as the unique measure satisfying 

(5 1 ,(J-z)- 1 5 1 ) = [ ^ zeC\R 
J R x — z 

where (•, •) denotes the inner product in £ 2 . It follows from the work of 
Janas and Naboko [8] that J\ ltVl has absolutely continuous spectrum 
covering the whole real line. 

l 
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This paper deals with random perturbations of Jai,?ji that are weak, 
in the sense that the variance of the perturbing random parameters 
grows like ~ n m with i] 2 <r]i- 

The case a(n) = 1, 772 < ("r^i = 0") and no perturbation 
off the diagonal, is the extensively studied family of one-dimensional 
discrete Schrodinger operators with a random decaying potential 
[3, 4, 5, 13, 14, 19]. For these Schrodinger operators it has been estab- 
lished that when r] 2 < — \, the absolutely continuous spectrum of the 
Laplacian is a.s. preserved. When i] 2 > — \, however, the disorder wins 
over and the spectrum is pure point with eigenf unctions that decay at a 
super-polynomial rate. At the critical point (7/2 = — \) the spectral be- 
havior exhibits a sensitive dependence on the coupling constant: The 
generalized eigenfunctions decay polynomially and the spectral mea- 
sure is pure point or singular continuous according to whether these 
generalized eigenfunctions are in £ 2 or not (for a comprehensive treat- 
ment of discrete Schrodinger operators with random decaying poten- 
tials, see Section 8 of [13]). 

From this perspective, the extension presented in this paper is in 
allowing growth of the off-diagonal terms. Intuitively, these terms are 
responsible for transport, and thus, their growth should have an effect 
on the spectrum similar to that of the decay of the potential (diagonal 
terms). Indeed, a particular case of our analysis is that of i]2 = 0, 
namely, the diagonal terms are i.i.d. random variables. We show that 
the critical point here is rji — |. Below this value, the spectrum is a.s. 
pure point, whereas above it the spectral measure is one-dimensional. 

More generally, let {X ul {n)}'^' =1 be a sequence of i.i.d. random vari- 
ables. Let {Yu){n)}°^ =l be another such sequence (the distributions of 
the X's and the Vs need not be the same). Assume the following is 
satisfied: 

(i) For all n 

(X u (n)) = (Y u (n)) = (1.2) 
(where (f(co)) = f n f{uj)dp and (Q, JF, dp) is the underlying probability 
space.) 

(ii) For any k G N 

(|r»| fe ) <oo, (|X»| fc )<oo. (1.3) 

(iii) 

<(X») 2 ) = 1, ({Y u) (n)f) = \ (1.4) 

(iv) The common distribution of X w (n) is absolutely continuous with 
respect to Lebesgue measure. 
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Given the quadruple T = (771, 772, Ai, A2) with < rji < 1, 772 < Vi 
and Ai, A 2 > 0, define the two random sequences 

br,U n ) = h M,mA n ) = A 2 n r?2 X aJ (n), a X2 , V2;u ,(n) = \ 2 n V2 Y UJ (n). 

Let 

ay,u>( n ) = a M,vi( n ) + a \2,m-A n ) 

and define 

■h,u> = J {{ar,U n )}n=i > { b rA n )}n=i) ■ ( L5 ) 
The assumptions on the parameters defining Jt,^, do not exclude the 
possibility that some of the off-diagonal terms will vanish. However, 
with probability one, this may happen only a finite number of times, so 
that Jx i£J has an infinite part with strictly positive off-diagonal entries. 
In the following, when we refer to Jr,u>, we refer to this part. 
We shall prove 

Theorem 1.1. For the model above, let 7 = rji — rj 2 an d let A = 
\ (A2/A1) 2 . The following holds with probability one: 

(1) 1/7 > I the spectrum ofJr >u > an d Hr^, the spectral measure 
of Jr,oj, is one- dimensional, meaning that it does not give weight 
to sets of Hausdorff dimension less than 1. 

(2) In the case 7 = | the spectrum of Jx,w is K and we have the 
following two possibilities: 

(a) If A > 1 — rji, then /i T w is pure point with eigenfunctions 
decaying like 

\^{n)\ 2 ~n-^\ 

(b) If A < I — rji, then /ix,w is purely singular continuous with 
exact Hausdorff dimension equal to 1 — j^z^ ■ 

(3) If 7 < I then the spectrum is pure point with eigenfunctions 
decaying like 

mn)\*~e-^- 2 \ 
In this case, if rji > 2r/ 2 then the spectrum fills R. 

Remark. We say that a measure, //, has exact Hausdorff dimension, g, 
if it is supported on a set of Hausdorff dimension g and does not give 
weight to sets of Hausdorff dimension less than g. For more information 
concerning the decomposition of general measures with respect to their 
Hausdorff-dimensional properties, consult [15] and references therein. 

Remark. In analogy to the Schrodinger case, one would expect to have 
absolutely continuous spectrum for 7 > 1/2. Unfortunately, though we 
believe this is true, one-dimensional spectral measure is all we could 
get. 
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Remark. The requirement that {X w (n)}^l 1 and {Y U! (n)}'^' =1 be iden- 
tically distributed sequences is not really necessary and is made here 
only to simplify the discussion. 

In resemblance of the Schrodinger case, the proof of this theorem 
follows by analyzing the asymptotics of solutions to the formal eigen- 
function equation " Jijj = Eip" . Namely, we shall analyze the solutions 
to the difference equation: 

a(n)ifj(n + 1) + b(n)ip(n) + a(n — l)ip(n — 1) = Eip(n) n > 1. (1.6) 

By a theorem of Kiselev and Last [12, Theorem 1.2], the results ob- 
tained have implications for the quantum dynamics associated with J; 
namely the behavior of a given vector tp under the operation of the one 
parameter unitary group e~ ltJ . More precisely, let X be the position 
operator defined by 

Theorem 1.2 of [12], which can be seen to hold in our setting, says 
that Theorem 1.1 and Proposition 2.7 below lead to 

Theorem 1.2. Assume 7 = | and A < 1 — 771 (namely, case 2(b) of 
Theorem 1.1). Then for any £ > 0, m G I, T > and ip G £ 2 , 

^ J J (x m e- itJr >"ij, e- itJr '-^j I dt > C(u, V>, m, e )T^ m /^-^' £ 

(1.7) 

with probability one. 

Note that r\\ may be chosen arbitrarily close to 1 while the spectral 
measure may have any dimension in [0, 1). Thus, by tuning the param- 
eters we obtain operators with any local spectral dimensions having 
arbitrarily fast transport. 

As an application of our general analysis, we study the Gaussian (5 
ensembles arising naturally in the context of Random Matrix Theory: 
The eigenvalue distribution functions for the three classical Gaussian 
ensembles are given by 

f PfN (E 1 ,..-,E N ) = J-e W (~J2 E j) II \ E i- E *f ( L8 ) 

l3N \ j=l J l<j<k<N 

with /3 — 1, 2 and 4 for the Gaussian Orthogonal Ensemble, Gaussian 
Unitary Ensemble and Gaussian Symplectic Ensemble respectively. 

A family of random matrix ensembles, indexed by j3, having fp^ as 
their eigenvalue distribution function, for any positive value of f3, was 
recently constructed by Dumitriu and Edelman [6]. The matrices in 
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these ensembles are finite random Jacobi matrices with the distribution 
of the off diagonal terms depending on (3: 

Definition 1.3. Fix (3 > 0. The random family of Jacobi matrices 
Jpv = J {{apA n )}n=i > ih^( n )}n=i) is defined by: 

(1) The random variables {ctp,u(n)}™=i , {b/3 t u(ri)}™ =1 are all inde- 
pendent. 

(2) bp }UJ {n) are all standard Gaussian variables (that is, with zero 
mean and variance= 1), irrespective of (3 and n. 

(3) The probability distribution function of ap^{n) is given by 

P{u | a^{n) <C} = -A- [ x^e'^dx. (1.9) 

1 {—) Jo 

In [6], Dumitriu and Edelman showed that the eigenvalue distribu- 
tion function of the finite matrix, obtained as the restriction of Jg iU; to 
the N x N upper left corner, is /^jv for any (3 > 0. 

From property 3 above, it follows that 



3 

n,2 



;i.io) 



{{aUn) ~ Mn))) 2 ) = \ + O . (1.11) 

Thus we see that the family J^ jW corresponds to the case rji = 1/2, 
r/2 = of the general matrices introduced above. Technically, the 
following theorem is not a corollary of Theorem 1.1, because of the 
O (n- l/2 ) term in (1.10) and the O (n^ 1 ) term in (1.11). The proof of 
Theorem 1.1, however, is robust with respect to such a change, and we 
have 

Theorem 1.4. For any (3, the spectrum of Jp iU) is R with probability 
one. 

If (3 < 2, then, with probability one, the spectral measure, /ip^, cor- 
responding to Jp^ and Si, is pure point with eigenf unctions decaying 
as 

\^u)(n)\ 2 ~ 

If ft > 2, then with probability one, for any e > 0, /j,p jUI has exact 
dimension 1 — | with probability one. Furthermore, for (3 > 2, we have 
that, almost surely, 



X m e~ itJ ^^, e- itJ e>"ij?) dt > C{u, 4>, m, e)T 2m ~ £ (1.12) 
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for any ip, e > and m. 

This result, without the dynamical part, was announced in [2]. We 
note that the analogous Circular j3 Ensembles can be realized as eigen- 
values of truncated CMV matrices. This was shown by Killip and 
Nenciu [10] and later used by Killip and Stoiciu in their analysis of 
level statistics for ensembles of random CMV matrices [11]. The bulk 
spectral properties of the appropriate matrices were analyzed by Simon 
[21, Section 12.7]. 

The proof of Theorem 1.1 is given in the next section. Since the 
proof of the spectral part of Theorem 1.4 is precisely the same, it is not 
given separately. As noted earlier, the dynamical part of our analysis 
(Theorem 1.2 and the corresponding statement in Theorem 1.4) follows 
immediately from Theorem 1.1 and Proposition 2.7, by Theorem 1.2 
of [12]. 

The method we use is a variation on the one used by Kiselev-Last- 
Simon [13, Section 8] in their analysis of the Schrodinger case described 
above. A notable difference is the fact that, due to the growth of 
the a(n), the effective energy parameter, ^|y, vanishes in the limit. 
This, in addition to requiring a modification in the technique of proof 
(see Lemma 2.4 and Proposition 2.5 below), leads to the fact that the 
asymptotics of the generalized eigenfunctions are constant over KL At 
the critical point ((771 — 772) = §), this implies uniformity of the local 
Hausdorff dimensions of the spectral measure. 

A modified Combes-Thomas estimate, for operators with unbounded 
off-diagonal terms, enters our analysis in the identification of the spec- 
trum of Jt,w Such an estimate may be of independent interest and 
thus is presented in the Appendix. 

2. Proof of Theorem 1.1 

We begin with a simple lemma that shows that, in a certain sense, 
Jt, w is a random relatively decaying perturbation of J\ um - 

Lemma 2.1. For any e > there exists, with probability one, a con- 
stant C = C{u),e) for which 

\n^X u (n)\ 



and 

K 2 F w (n)| 



< — (2.1) 



< — (2-2) 



where 7 = r\\ — r] 2 . 
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Proof. By (1.3) and Chebyshev's inequality we have for any k e N 

Pn = V | W 

By choosing 2/c > e _1 we see that 

oo 

J^P n < OO. 

n=l 

(2.1) follows now from Borel-Cantelli. The proof of (2.2) is the same. 

□ 

As stated in the Introduction, we follow the strategy of [13]. In par- 
ticular, we will deduce the spectral properties of Jt, w from the asymp- 
totics of the solutions to the corresponding eigenfunction equation. 

In order to fix notation, for a given Jacobi matrix 
J ({a(n)}™ =1 , {b( n )}™=i) an d a fixed E e R, denote by ip E a 
solution to the equation 

a(n)^ E (n + 1) + b(n)i/j E (n) + a(n - l)ip E (n - 1) = Eip E (n) n > 1. 

(2.3) 

It is customary to extend this equation to n = 1 by defining a(0) = 1. 
Clearly, the space of sequences {V' E ( n )}^o solving (2.3) is a two- 
dimensional vector space and any such sequence is completely deter- 
mined by its values at and 1. We let ip E (n) stand for the solution of 
(2.3) satisfying 

^J(O) = sin(0) = cos(0). (2.4) 

We note that, formally, ip E (n) satisfies 

J^ E = Eip E . 

Let 

/ E-b(n) q(n-l) 

S E (n) = »W "W 

and 

T £ (n) = S s (n) • - 1) • • • S E (1). 

Then, for any 0, 



and so 




\n^X u {n)\ 



> 



< 



C(k) 



n- 
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We call the matrices S E (n) defined above one-step transfer matrices, 
and for the matrices T E (n), we use the name n-step transfer matrices. 
Our main technical result is 

Theorem 2.2. Let J T u) be the family of random Jacobi matrices de- 
scribed in the Introduction. Then, for any E e R, the following holds 
with probability one: 



;i) ifi>\ 



(2) If 1 = 1 



(3) Ifl < \ 



y log II T*{n) || 2 (0 
hm \ — T~\ = 2 - 5 ) 



log II T E (n) II 2 
l im 1Qg tUL W I' = A - 771. (2.6) 



n^oo 



log(n) 



l im w& 11 " v ; 11 = A. (2.7) 

-,1-27 V > 



log llTf (n) 
n 1 

The EFGP transform (see [13]) is a useful tool for studying the as- 
ymptotic behavior of || T E (n) || in the Schrodinger case (a(n) = 1). For 
a(n) — > oo, certain modifications are needed. We proceed to present a 
version that is suitable for our purposes. 

Let J ({aui( n )}™=i{bui( n )}™=i) be a Jacobi matrix whose entries are 
all independent random variables. Let a{n) = (a u (n)) and = 
duj{n) — a(n), and assume that 

lim a{n) = oo (2.8) 

n— »oo 

and that 

hm ^# = (2.9) 
n^oo a(n) 

with probability one. These properties clearly hold for Jf tW (see Lemma 
2.1). In the analysis that follows we keep E e K fixed so we omit it 
from the notation. Define 



K u (n) 

Then, 



1 

a u (n) 



/ E-b u {n) 1 \ 

S u (n) = K u (n)S(n)K u (n - = a -( n \ a - (n) 

V M n ) J 

and 

f u {n) = S w (n) ■ S w (n - 1) • • • S u (l) = K u (n)T u (n). 
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For any 0, define the sequences {u (JJ ^{n)}'^ =1 and {vwA n )YnLi by 

twM \ = ft n \( cos O) 

v Ui<t> (n) J uK ' \ sin(0) 

so that, from the definition of Vty f° r t ne random Jacobi param- 

eters), we see that 

UuA n ) = Vw( n + !) *w( n ) = OwWV'w^W. (2.10) 

By (2.8) we see that for any E Gi and sufficiently large n, we may 
define k n e (0, 7r) by 

2cos(A; n ) = ^ T . (2.11) 

Clearly, A; n — > | as n — > oo. 

Now, define R^A 71 ) and O^A 71 ) through 

R^A 71 ) sin(0 w> 0(ra)) = u w ^(n) sin(/c n ) (2.12) 

and 

Ru>A n ) cos (9ojA n )) = a( n ) u uA n ) ~ v ojA n ) cos (kn) (2-13) 
so that (using (2.10)) 

RwA n ) 2 = v wA n ) 2 + a{ n ) 2u ioA n ) 2 ~ 2a(n)uwA n ) v "A n ) cos (^«) 
= ^{nf^A 71 ) 2 + a(n) 2 ^,0(n + l) 2 

which leads to 

Ru;A n ) 2 l _! 2a a; (n)a(n)^^(n) 2 

a(n) 2 (i,^(n) 2 + Vw( n + !) 2 ) «( n ) 2 (^W(n) 2 + tp u A n + J-) 2 ) 

| 2a^(n) 2 ^, (n) 2 
a(n) 2 (^,0(n) 2 + Vw( n + l) 2 ) 

Ed^n^A^^A 71 + 1) 
a(n) 2 (^^(n) 2 + ^^(n + l) 2 ) ' 

(2.14) 

By (2.9), it follows that the right hand side converges to one with 
probability 1, uniformly in 0, so that almost surely, for sufficiently 
large n, there are constants C 1: C 2 > such that 

CiR^in) 2 < a(n) 2 (^(n) 2 + ^ :<p (n + l) 2 ) < C 2 R^( n ) 2 

Now, by a straightforward adaptation of Lemma 2.2 of [13] it follows 
that for any two angles 0i ^ 2) there are constants C 3 , C 4 > such 
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that 

2 j ) /^\2\ ^ ;x/~,\2 II rri \ m2 



C 3 max(i2 u , ) 1 (n) ,R Wj( j >2 (n) ) < a(n) || T L 



[n\ 



2 2 ( 2 - 15 ) 

< C , 4 max(i? aJi 1 (n) , it^ 2 (n) ). 

Thus, we are led to examine the asymptotic properties of log R^A 71 )- 
Let us formulate a recursion relation for R^A 71 ) 2 '- (2.12) and (2.13) 
mean 

R^A 71 ) sin^^n)) 
R^A n ) coa(9 u ^(n)) 

sin(/c„) \ / u^A 71 ) 
a(n) -cos(k n ) J \ v w ^{n) 

sin(fc ri ) W 1 W ^{n + 1) 
a(n) -cos(k n ) J I a w (n) i I ^(n) 



We also know 



+ 2) \ _ , . / ^A 71 + 1) 



+ ; + I ^( n ) 



so 



-Ra;,0(n + 1) sin(6^(n + 1)) 
RuA n + l ) cos(^(n + 1)) 



sin(A; n+1 ) \ / 1 \ , , 

a(ra + l) -cos(A; n+1 ) J\0 a w (ra + l) ^ + i; 

1 \ 1 / sin(/c„) \ 1 / Ru,<t>{ n ) sin (^A n )) 



a w (n) y y a(n) - cos(/c n ) y \ Ru,A n ) cos (Q^A n )) 

, ° sin( ^ +i W>+i) 

a(n + l) -cos(fc n+ i) y v y 



sin(fc n ) \ 1 / R^A 71 ) sin %A n )) 



a(n) -cos(A; n ) y ^ R^A 71 ) C0S %A n )) 
Now, write 

g-ftq^fn+l) _ l 



(2.16) 



SJn+ 1) = a ? (n+1) x Mn+1) 
1 a w (n + l) 



£ i 



a(n+l) a(n+l) 



a(n + 1) 

a a; (n + 1) \ l a(n + 1) 

b^(n+l) q 



, a(n+l) 
"'" ^ Mn+l) 2 -a(n+l) 2 n 
a(ra+l) U 
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We define 

Sin(/c n+ i) \ / -^jy ~a(n+l) \ 

^ n + i >-[ a(n + 1) - cos(A;„ +1 ) ) { a(n + 1) J 

( sin(fc n ) VV sin(0 w> 0(n)) \ 
y a(n) -cos(A; n ) J \ cos(0 w> ^(n)) J ' 



and 



W(n + 1)-f sinC^+i) \( "%Sny 

KW + L) ~ \ a(n + 1) -cos(A; n+1 ) J \ a ^ n+ ^^ +1) 

_ / sin(£; n ) \ 1 / sin (6^,0(77,)) \ 
y a{n) -cos{k n ) J \ cos{O w ^{n)) J 

(we ignore the dependence on <fi since we keep it fixed). Then, from 
(2.16) we see that 

= i&Ltyi n Zu(n + 1) + w . (n + 1) r . (2 . 17 ) 

Oui,(f>(n) satisfies a recurrence relation as well: From (2.10) we have 
that 

VcuA n + 1) = M n + i)?w(n) (2-18) 

and 

0^(72 + 1)1x^(72 + 1) + 0^(72 + 1)7^,0(72) + a u (n)u Ui4> (n - 1) = Eu^n). 

(2.19) 

Write, using (2.12)-(2.13) 



COt(6^ 0^72 + 1)) = — r — 

sm{kn +1 )v u ^{n + 1) 
_ a(n + 1)7/^,0(72 + 1) - a^(n + 1) cos(A: ra+1 )M a;i0 (n) 
sin(A;„ + i)a w (n + 1)7^,0(72) 

(2.20) 

Furthermore, observing that 

Rw,<t>( n ) sin (^A n ) + M = 0(72) Sm(k n )u u ^(n), 

and 

RljA 71 ) cos(9 u ^(n) + k n ) = a(n) cos(k n )u u ^(n) - v^n) 

= a(n) cos(k n )u Ui4> (n) - 0^,(72)7^,0(72 - 1), 

we may write 

/ T1 \, l .\ aH cos(/c n )7/^0(n) - 0^(72)7/^,0(72 - 1) 

COt (6/^0 (72 J + fc„J = , . 1/./1J 

a(n) sm(A;„)^ i 0(72) 
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Substituting u L0 ,<j>{n + 1) from (2.19) into (2.20) and then u Uj< j,(n— 1) 
from (2.21) into the resulting equation, we get 

a(n + l)a(n) sin(A; n ) ,„ . . , . 
catflWn + 1)) = ^^-^ cat(Mn) + *») 

+ cot(k n+1 ) ( ~ a( ? + 1) l -l) 

V«-(n + l) 2 J (2.22) 

-bjn + l) 



sm(k n )a u (n + l) 2 
= /^(n + 1) cot(9 Ut4> (n)) + C(« + 1) 

where 

By (2.15) and picking two angles <pi ^ 02, Theorem 2.2 follows from 

Proposition 2.3. Let Jy t0J be the family of random Jacobi matrices 
described in the Introduction. Then, for any £gR, and for any (ft, the 
following holds with probability one: 

(1) Ifl>\ 

log RE An) 2 
lim ° ; = m (2.23) 

n->oo log(n) 



(2) //7=| 



(3) //7 < \ 



log R E An) 2 

k4rr- = AH- (2.24) 

n^oo log(n) 



logiff» 2 



hm ° = A. (2.25) 



Proof. As in [13] we shall prove the statement by using the recursion 
relation for Ruin) 2 (equation (2.17)). Namely, we shall prove that 



± (log (|| ZJj) + WJj) f) - log ) (2.26) 

converges to the appropriate limit, where F^(n) = log(n) for 7 > | 
and Fj(n) = n 1-27 otherwise. From this point on, a u (n) = ay t u(n) and 
b u (ri) = b r ,U n )- 

We shall need some estimate on the behavior of 9^{n). We start with 

Lemma 2.4. For any e > 0, there exists, with probability one, a con- 
stant C = C{uj,s), such that 

\e Ut4> {n + l)-e Ut4> {n)\ < Cmax^.n" 1 ). (2.27) 
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Proof of the Lemma. We start by proving a similar statement for 
Mn + l)-l| + |C(n+l)|: 

Recall that k n — > |. Thus, for large enough n, sin(/c n ) > \. More- 
over, cos(/c„) ~ It follows that 



| sin(/c„) - sin(A;„ + i)| 

Thus, 

|^(n+l)-l| + |C(n + l)| < 



cos(A; n ) 2 - cos(/c n+ i) 



< 



C" 



sin(fc n ) + sin(A; n+1 )| ~ n l + 2 ^ 
sin(/c„) / a(n + l)a(n) 



+ 
+ 
+ 



sin(/c n+ i) V a w (" + l) 2 
sm(k n ) - sm(k n+1 



cot (/a. 



sin(/c„+i) 

a(n + 1) 



n+lj 



a^in + l) 2 
a(n + 1) 



sin(/c n )a^(n + l) 2 
= h + h + h + h- 
Now, choosing e < r\\ — r/ 2 , by Lemma 2.1 



6 w (n + l) 



7i < 2 



a(n + 1) (a(n) — a(n + 1)) 



+ 2 



a w (n + l) 2 
a w (n + l) 2 



4 



ct^in + l)a(n + 1) 



h < 2 1 cos(/c n+ i 



a w (n + l) 2 
< Ci(lu,e) max(r?T 7+e , n^ 1 ) 



a u {n + l) 2 

2a u] {n + l)a(n + 1) + a u (n + l) 2 



a u {n + l) 2 



< 



and 



7 4 < 2 



a(n + l)b u (n + 1) 



a w (n + l) 2 

almost surely. We have shown above that 

C 2 



< 



h < 



n 



so we see that there exists, with probability one, a constant C 5 (cu,e) 
for which 



Mn + 1) - 1| + |C> + 1)] < -^(n" 7+£ , r^ 1 ). 

max 



(2.28) 



Now, use 



1 - 



l + i cot a; 
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and (2.22) to see that, if ^(n + 1) - 1| + 2|C«( n + 1)| < |, (which 
indeed happens almost surely, for large enough n), then 



<4(Mn+l)-l| + |C w (n+l)|), 



which implies (by \e lx — 1| > — ) that 

KA n + !) - ^Wl < ^ (M™ + !) - !| + IO + 1)1) ■ 

This, together with (2.28), implies (2.27) and concludes the proof of 
the lemma. □ 

The direct consequence of this is 
Proposition 2.5. Assume f(n) is a function that satisfies 

f(n) = o(n r ) 



and 



with 



Then, for any (f>, 



f(n + l)-f(n)=o(n r - 1 ), 



7-1 if 7 > 1/2 and r] 2 > 

—7 if 7 < 1/2 and r\ 2 > 

7)i — l if 7 > 1/2 and i] 2 < 

7) 1 - 27 if 7 < 1/2 and r/ 2 < 0. 



lim p7=tE/0') oo8(2^(j')) = 



(2.29) 



(2.30) 



almost surely. The same statement holds with 9 replaced by 9. 



Proof of the Proposition. We shall prove the statement for 9. By sum- 
mation by parts, 



£/O')cos(20 w ,,O')) 



3=1 



n-l j 



/(n) E cos(20 w> ,(j)) "EE «*(2*«,*(0) (/(J + 1) " /(J)) • 
j=i i=i «=i 

(2.31) 
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Thus we are led to examine Yla=i cos (2#c,</>(0)- Assume that j = 2m 
is even. Then 



5>os(2^(Z)) 



1=1 



= (cos(20 Wi0 (2Z)) - cos(20^(2Z - 1) + tt)) 
z=i 

< l 2 ^(20 - 2^(2/ — 1) — tt| 



z=i 



< 2^ |^(2/)- ^(2/ -r 



7T 



"2Z-1 



< max(j 



almost surely, by Lemma 2.4 and by 



1— 7+e -1—T/i 



(2.32) 



7T 



< 2 1 cos(fc n )|, 



which holds for sufficiently large n. Thus, for any j, we get 



^cos(2^(0) 



< ^max^-^j- 1 -" 1 ) + 1 
<(C w + l)max(j 1 - 7+e ,j 1 - ,n ). 



(2.33) 



A simple calculation finishes the proof for 0. The proof for 9 follows 
the same argument, with an additional \k 2 i — foj-i] term in (2.32). □ 



We abbreviate 



A ^ = Mi) 2 - a(j) 2 = 2 Mj) + ojjf 



a(j) 2 



Ai 



Ai 



(2.34) 



J 



2-> 



and 



S(j) Ai p 
By a straightforward calculation we get 

|2 sin 2 (^(j - l))S(j) 2 



(2.35) 



sin 2 (%_!)a(j - l) 2 

x I A„(j) 2 + B^jf + 2cos{k j )B UJ {j)A lJJ {j) 



(2.36) 
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II ^0) || = 1 + sm (^(j - 1))- 



+ 



~a(j - l) 2 

5(j) 2 sin 2 (^(j - 1)) sin 2 (A;,) - sin 2 (Vi) 



(2.37) 



a(j - l) 2 sin 2 (%_i) 
and 



5(j - 1) 

x l ^Bin(^ -i)) / 

sm 2 (A; J _i) Va(j-l) 



- cos(fcj) sin(0 w (j - 1))) - 



^(j>in(2^(j - 1) ) 
2 sin(A; J _ 1 ) 

(2.38) 



(Recall 8 u (j) = 6^{j) = 9^{j) + kj = u {j) + k j .) 

Let e < |. By Theorem 2.1, the fact that cos(kj) ~ and the 
identity 

| sin 2 (a;) — sin 2 (y)| = | cos 2 (re) — cos 2 ' 



it follows that || W u (j) H 2 = a(j~ 27+2£ ), {ZM,WM) = O w {j-^) 
and (|| Zu{j) || 2 —1) = 0(j~ l ) with probability 1, where the notation 
O u indicates that the implicit constant depends on uj. Thus, we can 
use 

x 2 

log(l + x) =x- — + 0(x 3 ), 
together with the observation that (almost surely) 

( || zm) II 2 -i + 2(^(j),w.(j))+ II WM II 2 ) 2 

= 4(^(j),W.(j)) 2 + a(-^z7 
to see that, with probability one, for large enough j we have that 

log (l+ || Z w (j) || 2 -1 + 2(Z 0J (j), W u {j))+ II w w (j) II 2 ) 

= (|| ZUj) || 2 -1) + 2(ZUj), W w (j'))+ II W w (j) || 2 (2.39) 
-2(WW w (j')) 2 +a[^: + 



i'37-e 9 'l+7-s 
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Therefore, since the edition of a finite number of terms is inconsequen- 
tial, it follows that 

i - 



£«^)£Mll^) + nu,)f) 



3 

n 



lim pT^tE (dl Z ^ 'I' -1) + 2(^(j), II W w (j) 



J'=l 

2(3,0'), W w (j)) 2 ) 



(2.40) 



with probability one, in the sense that both limits exist together and 
are equal if they do. 
Similarly, 

J™ E ^ (|f ) - >- j^j E "og (i + 

™ F 7 (n) ^ V 2 / 

(2.41) 

Thus, our problem is reduced to computing the limits: 

n 

& s "»j^E(ll*0>ll'-i) 

^ V ' 7=1 



1 n 

£w = lim-— -^HW^i) || 2 

tzy,^ hm ^ ^ 
—2 " 
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By (2.37), 



ZM II -1 = sin (9 u {j - 1))- 



+ 



a(j - 1) 2 

a(Jf su^jejj - 1)) sin 2 (%) - sin 2 ( y x ) 
a(j - l) 2 sin 2 



The last term on the right is absolutely summable (= G(n 1 2r/1 )) so 
we only need to look at sin 2 (^(j - l)) 5 ^^ 1 ^ . But, by sin 2 (a) = 

\ — cos ^ 2a - > and by Proposition 2.5, we see that (recall a(j) = \ij~ vi ) 

, _1 1 a(j) 2 ~ ajj ~ I) 2 

^ 2 n^oo F 7 (n) 4^ «(J - I) 2 

J (2.42) 

- 1 hm 1 V 2Vl - f * if ^^V2 

Z n— >oo 

F 7 (n) ■f-' j \ otherwise. 



For the other four limits, we shall use extensively Lemma 8.4 of [13], 
in order to replace A u and B w by their means. For example, write 



(2.43) 



Then the first three terms have mean zero and so, by Lemma 8.4 of 
[13] we get that 
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Expanding (A^j) 2 ), throwing out terms that are o(j 27 ), and applying 
Proposition 2.5 (writing, again, sin 2 (o;) = | — cos ^ 2a ) ) ; we g e t 



1 >r-^ j 2m 



W^j-l) 2 ^ 2 ^)' 
x «4YL(j) 2 > + <X w (j) 2 » (2.44) 

= lim^V2Ar 2 ^ = l ° if7>V2 
n^oo FJn) ^ I 2A otherwise, 

rv j =1 v. 

(recall (K,(j) 2 ) = |). Applying the same procedure to £ Z w an d £a we 
get 



f , 2 " g» sin 2 (^(j-l)) / oa F4j) 2 

^ ~ *>) ^ (j -l)-n sin 2 (Vi) \ J 2 ^ 

= lim 1 V A 1 - f ° if 7>V2 
n^oo F 7 (n) 4^ 2 j 2 t \ f otherwise 



(2.45) 



and 



n^oo F» ^ \ J 27 



1 A/A 1\_ f if 7 >l/2 



™ F 7 (n) V 2 f y J If otherwise. 



(2.46) 



The computation of £zw 2 involves sin 4 {9 J) — | — cos ( 26> ") _|_ cos (4^) ^ f Qr 
which Proposition 2.5 is useless. Luckily, the cos{A9 w {j)) cancels out. 
As before 



x 



+ 



( 4^(j) 2 )sin 4 (^(j-l)) 

sin^^-i) (j - l) 2 " 1 

(X w (j) 2 )sin 2 (2^(j-1))~ 
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Now write 

'(4^(j ) 2 )sin 4 (^(j-l)) {XUm^y / I )) 

sin 2 ^) (j - l) 2 ^ 

's in 4 (fl„(j-l)) j 2 ^ sm 2 (2^Q--l) )" 
sin 2 ^-^) (j - l) 2 ^ + 4 



8(j - l) 2 ^ sin 2 (%_ 1 ) 8 

? -2r?i 

+ 20--.)l m -( tH ) mW - 1|) 

+ cos(4^ -l))( ( ._ 1)2 f^ fe _ ir l 
3?' 2t?1 1 

+ 



8(j-l) 2 ^sin 2 (%_ 1 ) 8 

+ or 1V 4 • 277 ; cos(20 w (j - 1)) 
2(j - l)*n sin 

+ cos(40 w (j " !)) ^fa -1 ) + 0(^ 2m ) 
to see that 

, _f if 7 > 1/2 , , 

kzm - | _ 2A otherwise; 

where the cos^^) term vanishes by Proposition 2.5. Summing up the 
various limits, the proposition is proved. □ 

Proof of Theorem 2.2. By (2.15), the theorem follows from Proposition 
2.3. □ 

Theorem 1.1 almost follows immediately from Theorem 2.2 and 
Proposition 2.3. As in the Schrodinger situation, the case 7 = \ re- 
quires some subtle reasoning. We have established that, in this case, 
with probability one, equation (2.3) has a solution, ip, with 

l^njpxn^ 1 . (2.48) 

In order to use subordinacy theory, we need the existence of another 
solution with faster decay at infinity. The following is Lemma 8.7 of 
[13], formulated for general regular matrices: 
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Lemma 2.6. Let u$ = (cos 0, sin 0) G M 2 . For any matrix, A G 
GL 2 (M) with det(A) = d > 0, let (p(A) be the unique (p G (-§, f ] 
^|| ^ ||=V5||A||-i. Define p(A) = M^L. 

Let A n be a sequence of matrices in GL 2 (M) with det(A n ) = d n > 0, 
that satisfy -k= II A n II— > oo and cL f f "n, 1 ^" I — > as n — > oo. Lei 

vUT 11 n 11 n ||A„||j|A n+ i|| 

p n = p(A n ) and (p n = (f)(A n ). Then: 

(1) (p n has a limit (poo if and only if 'Hindoo p n = p^ exists (p^ = oo 
is allowed, but then we only have \(f> n \ — > |). 

(2) Suppose (p n has a limit (p^ ^ 0, | (equivalently, p^ ^ 0,oo). 
Then 

log || A n Moo || -\ogy/a\ 

hm , = — 1 



n^oo log || A n || - log 

if and only if 

v log |p»- Pop) 

limsup -j= < —I. 

n->oc log || A n || -logVdn 

Proposition 2.7. Let Jt, w fre £/ie family of random Jacobi matrices 
described in Theorem 2.2, with 1 = \- Then, for any E G K, £/iere 

exists, with probability one, an initial condition VP^a,) = ^ gin(0(u;)) 

Um log II II = 1 

log n 2 



Proof. We imitate the proof of [13, Lemma 8.8]. By Proposition 2.3 

\Rw,o{n)\ = 1 
logn 2 



n^oo 

and 



lim log|^,o(n)| = l (A + ^ 



log | R^/2{n) | 1 

hm = - (A + 771) 

log n 2 



for almost every tu. By (2.12)-(2.13), 

Ru,,o( n ) R w,K/2(n) sin(6> Wi7r / 2 (n) - w ,oW) = fi H sin£;„ 
so that, for almost every cu, 

lim lQ gl^oW-WWI = _ A (2 5Q) 

n->oo log n 

(recall d(n) = \in m ). 

Let A-W = Tg^gy- Then 

L w (n) =logp a ,(n + l)-logp w (n) = log(l + A , Wi0 (n)) -log(l + A^ i7r/2 (n)) 

(2.51) 
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X Ut4 ,(n) = (|| Z u j(n + 1) + W Ut4> (n + 1) f -l) . (2.52) 

Since Xui,cf>(j) — O {-j^J almost surely, we may apply a finite Taylor 
expansion to the above (using (2. 36), (2. 37) and (2.38)) to see that, 
with probability one, for large enough n 

L u {n) = Al(n) (sin 2 (0 W)O (n - 1)) - sin 2 (^ i7r/2 (n - 1))) 

+ A 2 (n) (sin(2^ i0 (n - 1)) - sin(2^ )7r/2 (n - 1))) 

+ Al(n) ( fw(9 Ut0 (n - 1)) 8m% t0 (n - 1)) 
- sin(^ i7r/2 (n - 1)) sin(0 Wi7r / 2 (n - 1)) 



+ e>( 



-(l+A-e) 



) 



where 
and 



<Ai(n)> = (Al („)) = <A^ („)> = 



(a>)) 2 ) < 



c 



i = 1,2,3. (The O (n ( 1+A £ - ) ) for the remainder follows from Lemma 
2.1 and (2.50).) 

A standard application of Kolmogorov's inequality and the Borel 
Cantelli Lemma shows that with probability one, for large enough k, 
and for each i = 1,2,3 



sup 

m=2 k ~ 1 +l,...,2 k 



E A -c?) 

j=2 k - 1 +l 



< k 



and also 



sup 

m=2 fc - 1 +l,...,2*-l 



< fc. 



(2.53) 



(2.54) 



Combining this with the fact that, with probability one, for large 
enough n, 

(sin 2 (0 W)O (n - 1)) - sin 2 (^ i7r/2 (n - 1))) = O {n^) , 
(sin(2^, (n - 1)) - sin(2^ i7r/2 (n - 1))) = O (n- A+£ ) 

and 

sin(^ i0 (n - 1)) sin(6> Wi o(n - 1)) 
- sin(^ /2 (n - 1)) sin(^ /2 (n - 1))) = O (n- A+£ ) , 
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it follows that 



exists and 



n=l 



5> 

n=N 



n 



< C ul N- K+£ 



almost surely. 

Thus, with 
Puj (oo) exists and 



Thus, with proability one, lim n ^ 00 B fl "' o( ?\ = lim n ^ 0O p a; (n) = 

Il uj,ir/2\ n j 

v log\p u (n) -p w (oo)| 
hmsup < —A. 

n-oo log U 

Lemma 2.6 completes the proof (note that d n = det T^(n) = - j^ )- D 
We are now ready to complete the 

Proof of Theorem 1.1. (1) By Theorem 2.2, Fubini's Theorem, the 
fact that the distribution of X u (n) is absolutely continuous with 
respect to Lebesgue measure, and the theory of rank-one pertur- 
bations ([20]), it follows that, with probability one, the spectral 
measure is supported on the set of energies where, for any e > 
and sufficiently large n, || T^{n) || 2 < n~ m+£ . From Corollary 
4.4 of [9] it follows now that the spectral measure is contin- 
uous with respect to (1 — £)-dimensional Hausdorff measure, 
for any e > 0. Thus the spectral measure is one-dimensional. 
Since ip(n) = (5i, (J — z)^ 1 5 n ) solves the eigenvalue equation 
for z (away from n — 0), Theorem A.l and Wronskian conser- 
vation imply that if z G K were to be outside of the spectrum, 
the transfer matrices would have to exhibit exponential growth. 
Since this is not the case, it follows that the spectrum is R. 

(2) In this case again, the fact that the spectrum is R follows from 
the polynomial bound on the transfer matrices in Theorem 2.2 
and Theorem A.l below. As for the properties of the spec- 
tral measure, these follow from Theorem 1.2 in [9] using (2.48), 
Proposition 2.7 and the theory of rank one perturbations. 

(3) The existence, with probability one, of an exponentially decay- 
ing eigenf unction, for every E e R, follows from Theorem 2.2 
and Theorem 8.3 of [16]. Fubini and the theory of rank one per- 
turbations imply that the spectral measure is supported, with 
probability one, on the set where these eigenfunctions exist. 
Comparing powers of n in the exponent, Theorem A.l implies 
that, as long as r]i > 2r/ 2 , the spectrum fills R. 
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□ 

Appendix A. A Combes-Thomas Estimate for Jacobi 
Matrices with Unbounded Parameters 

This section presents a Combes-Thomas estimate, suitable for appli- 
cation to Jacobi matrices with unbounded off-diagonal terms. Also see 
[18] for a related result. 

Theorem A.l. Let J = J ({a(n)}^ =1 , {b(n) }™ =1 ) be a self-adjoint Ja- 
cobi matrix such that < a{n) < f(n) for a nondecreasing function 
f(n). Let z G C be such that dist(z, Spec( J)) = a (where Spec(J) is 
the spectrum of J). Then 

\(6i,(J-z)- 1 6 N )\<-e- a »- N (A.l) 

a 

where = min ^1, j^j^j ■ In particular 

| (Jr,. - z)' 1 5 N ) | < ^ e - c ^> Nl - m (A.2) 

almost surely. 

Remark. The monotonicity of / is not essential. One may instead 
replace f(N) in the formula for ckjv, by max(/(l), . . . , f(N)). 

Proof. Let Rn be the diagonal matrix defined by 

r> / \ f eQJV ' n n < N , . 

RN{n,n) = | eajv .jv n > ^ (A.3) 

Then 

e a„.(AT-l) (J _ ^ = ^-1 ( J _ ^-1 

so it suffices to bound ji R N l (J — z)' 1 Rn ||. Noting that 
i?^ 1 (J — z) 1 -Rat = (Rjj 1 (J — 2) Rn) = Cn(z), we may apply the 
Resolvent identity to get 

C N (z) = {J- z)- 1 + C N {z) -(J-z-R^(J- z) R N ) (J - z)- 1 . 

A simple computation shows that || J — z — R^ 1 (J — z) Rn \\< 
2ef{N)a N if a N < 1. Thus, by || (J - ||< i , we see that 

\\C N (z) \\< I + C N (z)^^-a N 
a a 

so, by aN < \ • 2ef(N) ' we see ^ na ^ 

II Cat(^) ||< - 

(7 
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which finishes the proof. □ 
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